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ABSTRACT 

We present a strong-lensing analysis of the galaxy cluster MACS J1206. 2-0847 (z=0.44) using UV, 
Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, 
with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, 
combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant 
sources. These images, along with the previously known arc, span the redshift range 1 < z < 5.5, and 
thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, 
the inner mass-profile of this cluster. We find an inner profile slope of d log S/d log 8 ~ —0.55 ± 0.1 
(in the range [1", 53"], or 5 < r < 300 kpc), as commonly found for relaxed and well-concentrated 
clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii 
for different source redshifts. For a source at z s ~ 2.5, the critical curve encloses a large area with an 
effective Einstein radius of 9e = 28 ± 3", and a projected mass of 1.34 ±0.15 x 10 14 M©. From the 
current understanding of structure formation in concordance cosmology, these values are relatively 
high for clusters at z ~ 0.5, so that detailed studies of the inner mass distribution of clusters such as 
MACS J1206. 2-0847 can provide stringent tests of the ACDM paradigm. 

Subject headings: dark matter, galaxies: clusters: individuals: MACS J1206. 2-0847, galaxies: clusters: 
general, galaxies: high-redshift, gravitational lensing: strong 



1. INTRODUCTION 

Massive galaxy clusters, due to their high inner mass- 
density, are known to form prominent gravitational 
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lenses. The expected distribution of lens-sizes and the 
abundance of giant lenses in particular, have been now 



established by N-body simulations (e.g., [Hcnnawi et al 
2007), sem i-analytic calculations (e.g., Oguri & Bland- 
ford 2009), and recently, also examined observationally 
by a statistical a nalysis and lens-m odeling of 10,000 
SDSS clusters (see |Zitrin et al.||2011b ). 

Due to the hierarchical growth of structure in the 
Universe, collapsed, virialized clusters should be found 
mostly at lower redshifts. These clusters make for ex- 
cellent lenses as there is more mass concentrated in the 
cluster center, boosting the critical lensing area. Ac- 
cording to this assumption, along with the dependency 
on the cosmological distances involved, lensing should be 
therefore optimized in clusters at redshifts of zi ~ 0.2. 
However, recent work has uncovered more large higher- 
redshift (zi ~ 0.5) lenses than expected by ACDM and 
relate d simulations, even after accounting for lensin g bias 



(e.g., Zitrin et al.||2011a| 
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history of growth, or a non-Gaussian distribution of mas- 
sive perturbations. 

The galaxy cluster MACS J1206.2-0847 0=0.4385; 
MACS1206 hereafter), is an X-ray selected system at in- 
termedi ate redshift found by the Massive Cluster Survey, 
MACS ( |Ebeling et aT1[200r] [20071 [20To| >, and therefore 
constitutes an interesting lensmg target. A first mass 
model for this cluster was presented by |Ebeling et al. 
(2009), based on 1-band HST/ACS imaging (F606W), 
combined with additional opti cal and NIR ground-based 
imaging. Ebeling et al. ( 2009[ ) have identified one multi- 
ple system, consisting of a giant arc and its counter image 
at z s = 1.036, and presented a mass distribution for this 
cluster, though without constraining the profile due to 
the lack of sufficient high-resolution color-imaging, and 
correspondingly, other multiple-systems. The 16 HST 
bands chosen for the the Cluster Lensing And Supernova 
surve y with Hubble (CLASH) project ( [Postman et al. 
2011| , ranging from the UV through the optical and to 
the IR, along with spectra from the VLT/VIMOS for 
some of the brighter arcs, enable us to obtain accurate 
redshifts for the multiply-lensed sources presented in this 
work. We use these data available to date, alo ng with 
our well-tested approach to SL modeling 



hurst et al.||2005[ |Zitrin et al.||2009a|b| l2010[ |20ira|cf 
to find a significant number of multiple images across the 
central field of MACS1206 so that its mass distribution 
and inner profile can be constrained for the first time, 
and with high precision. 

The paper is organized as follows: In §2 we describe 
the observations, and in §3 we detail the SL analysis. 
In §4 we report and discuss the results. Throughout we 
adopt a concordance ACDM cosmology with (fl m o = 0.3, 
^ao = 0.7, h = 0.7). With these parameters one arcsec- 
ond corresponds to a p hysical scale of 5.67 kpc for this 
cluster (at z = 0.4385; |Ebeling et al~1|2009| ). The refer- 
ence center of our analysis is fixed on the brightest clus- 
ter galaxy (BCG): RA = 12:06:12.15 Dec = -08:48:03.4 
(J2000.0). 

2. OBSERVATIONS AND REDSHIFTS 

As part of the CLASH program, MACS1206 was ob- 
served with HST from 2011 March to 2011 July This is 
the third of 25 clusters to be observed to a depth of 20 
HST orbits in 16 filters with the Wide Field Camera 3 
(WFC3) UVIS and IR cameras, and the Advanced Cam- 
era for Surveys (ACS) WFC. The images are processed 
for debias, flats, superflats, and darks, using standard 
techniques, and are then co-aligned and combined us- 
ing drizzle algorithms to a scale of 0.065"/ pixel. The 
full UVIS/ACS/WFC3-IR data set is then importantly 
used for multiple- images verification and measurement of 
their photomet ric redshifts using both the BPZ pro gram 
(|Benftez||2000] | Bemtez e faLJ[2004l |Coe et al. |[2006| ) , and 
LePhare ( |Arnouts et al.||l999| |llbert et al.||2006p , where 
in practice 15 bands were used for the photometry, as 
observations for this cluster were still in progress during 
the preparation of this paper (nearly all 20 orbits had 
been completed for most filters. The F3 36W band was 
not us ed) . Further details are presented in Postman et al 

(poll) 



spectroscopy of 14 southern CLASH clusters, targeting 
hundreds of cluster members per cluster and SL features 
in their cores. Details on this program will be presented 
elsewhere, when observations for one cluster are com- 
pleted. For each cluster, four VIMOS pointings are used, 
keeping one of the four quadrants constantly locked on 
the cluster core, thus allowing long exposures on the 
arcs, where exposure times for each pointing are about 
45-60 minutes. By filling the inter-quadrant gaps, the 
final VIMOS layout covers 20-25' across. Either the low- 
resolution LR-Blue grism or the intermediate resolution 
orange MR grism is used, depending on the photometric 
redshifts of the targets. 

The spectra presented here are the results of the very 
first observations for this program, consisting of four 
pointings with the LR-Blue grism, obtained in 2011 
March-April, which yielded approximately 1000 red- 
shifts. This configuration provides a spectral resolution 
of ~ 28A with 1" slits and a useful wavelength coverage 
of 3700-6800A. Preliminary HST/CLASH data from the 
first two visits of MACS1206 were used to select images 
1. [1,2,3], 2. [1,2, 3], 3. [1,2, 3], and 4.1, as spectroscopic tar- 
gets (see Figure fl]). The slits ran along the NS direction. 
For some exposures the seeing was very good so that sep- 
arate spectra of the pairs 2.1/3.1 and 2.2/3.2 (separation 
of ~ 1") were taken (and are shown in Figure [2] Table 
[I]), though in other cases the pairs were blendea in the 
slit. The spectrum of 1.2, which contains the blend of 
the giant arc and a compact cluster galaxy, is not shown 
A spectrum of this arc, covering redder wavele ngths in- 
cluding the [Oil] line, was however published in Ebeling 



We obtain spectra for a number of multiple systems 
uncovered here, taken as part of the VLT/VIMOS Large 
Programme 186.A-0798, which will perform panoramic 



et al. (2009). Our spectroscopy also confirms that the 
tour compact galaxies right on the East and West side to 
the giant arc are early-type cluster members. 

3. STRONG LENSING MODELING AND ANALYSIS 

The approach to lens modeling we us e here (e.g., 
Broadhurst et al. 2005 Zitrin et al. 2009b) begins with 
the assumption that mass approximately traces light, so 
that the photometry of the red-sequence cluster member 
galaxies is used as the starting point for our model. In 
particular, we use the F814W and F475W bands to filter 
in the brighter member galaxies (tosi4 < 23 AB mag), 
and the F814W flux to derive the relative weight of each 
member. Using the extensive multi-band imaging and 
corresponding photometric redshifts, these galaxies can 
be then verified as members lying at the cluster redshift. 

We approximate the large scale distribution of clus- 
ter mass by assigning a power-law mass profile to each 
galaxy, the sum of which is then smoothed, using a 2D 
spline interpolation. The polynomial degree of smooth- 
ing (S) and the index of the power-law (q) are the most 
important free parameters determining the mass profile: 
steeper power-laws and higher 2D polynomial de grees, 
generally entail a steeper profile (see Zitrin et al. 2009b). 
A worthwhile improvement in fitting the location of the 
lensed images is generally found by expanding to first 
order the gravitational potential of this smooth compo- 
nent, equivalent to a coherent shear describing the overall 
matter ellipticity. The direction of the shear and its am- 
plitude are free parameters, allowing for some flexibility 
in the relation between the distribution of dark matter 
(DM) and the distribution of galaxies, which cannot be 
expected to trace each other in detail. The total de- 



CLASH Mass Profile of MACS J1206.2-0847 
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Fig. 1.— Galaxy cluster MACS1206 (z = 0.4385) imaged with HST/ACS/WFC3. North is up, East is left. We number the multiply- 
lensed images used and uncovered in this work. The numbers indicate the 50 lensed images, 47 of which are uncovered here and correspond 
to (at least) 12 newly identified sources and candidates, and the different colors are used to distinguish between them. Note that candidate 
systems are marked in "c" . Details on the each system are given in Table [l] The overlaid white critical curve corresponds to system 4 at 
z s = 2.54, enclosing a critical area with an effective Einstein radius of ~ loTTkpc at the redshift of this cluster (28"). Also plotted is a blue 
critical curve, which corresponds to the lower redshift of system 1, the giant arc system at z s = 1.033. The composition of this color image 
is Red=F105W+F110W+F125W+F140W+F160W, Green=F606W+F625 W+F775W+F814W+F850LP, and Blue=F435W+F475W. The 
bottom-left inset shows an enlargement of the central core. 



flection field 3t(0), consists of the galaxy component, 
&gai(0), scaled by a factor K ga i, the cluster DM compo- 
nent o.dm{6), scaled by (1-K ga i), and the external shear 
component a ex (9): 

a T {0) = K gal a ga i{6) + (1 - K ga i)a DM {0) + a ex (9). (1) 

The best fit is assessed by the minimum \ 2 uncertainty 
in the image plane, and the errors are determined accord- 
ingly, by adopting a positional error of 2". This we have 
found is a typical value, following previous findings for 
the e ffect of large-scale structure (L SS) along the line-of- 
sight (Julio et al. 2010 Host||201lj ), and by propagating 
a typical Az ~ 0.1 photo- z uncertainty (for z s ~ 2) 
into the image-plane reproduction error. Note also, this 
image-plane minimization does not suffer from the bias 
involved with source-plane minimization, where solutions 
are biased by minimal scatter towards shallow mass pro- 
files with correspondingly higher magnifications. 



It should be stressed that the multiple-images found 
here are accurately reproduced by our model (e.g., Fig- 
ure [3]) , and are not simple identifications by eye. Due to 
the small number of parameters in our model, it is ini- 
tially well-constrained, enabling a reliable identification 
of other multiple-images in the field. The mass model 
predictions are identified in the data and verified further 
by comparing the SEDs and photometric redshifts of the 
candidate multiple-images. The model is successively re- 
fined as additional sets of multiple images are incorpo- 
rated to improve the fit, importantly using also their red- 
shift information for better constraining the mass slope 
through the cosmological relation of the Di s /D s growth. 

4. RESULTS AND CONCLUSIONS 

Using CLASH imaging, we have identified 47 new 
multiple-images in MACS1206, corresponding to 12 dis- 
tant source s. These images, a long with the previously 
known arc ( Ebeling et al.||20"09 ), span the redshift range 
1< z< 5.5 and thus allow us to derive a robust mass 
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TABLE 1 

Multiple-image systems and candidates 



ARC 


RA 


DEC 


BPZ z p h ot 


LePhare z phot 


spec-z 


Z-model 


A Position 


Comment 


ID 


(J2000.0) 


(J2000.0) 


(best) [95% C.L.] 


(best) [90% C.L.] 






(arcsecs) 




1.1 
1.2 
1.3 


12:06:10.75 
12:06:10.82 
12:06:11.29 


-08 

-08 
-08 


48:01.01 
48:08.95 
47:43.44 


-H 

1.04 [0.96-1.12] 
1.01 [0.93-1.09] 


-H 

1.06 [1.04-1.09] 
1.05 [1.03-1.06] 


1.033 

(1.036) 

1.033 


(1.033) 


3.3 
0.1 
3.3 
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1 2009 I 








2.1 


12:06:14.53 


-08 


48:32.37 


-H 


-H 


3.03 


(3.03) 


0.2 




2.2 


12:06:15.00 


-08 


48:17.67 


3.40 [3.23-3.57] 


3.19 [3.12-3.36] 


3.03 


" 


0.5 




2.3 


12:06:15.03 


-08 


47:48.07 


3.68 [3.50-3.86] 


3.64 [3.59-3.70] 


3.03 


" 


2.1 




3.1 


12:06:14.43 


-08 


48:34.20 


3.73 [3.55-3.92] 


3.65 [3.60-3.73] 


3.03 


(3.03) 


0.2 




3.2 


12:06:15.00 


-08 


48:16.50 


-H 


-H 


3.03 


" 


0.5 




3.3 


12:06:15.01 


-08 


47:48.65 


3.52 [3.34-3.70] 


3.62 [3.53-3.67] 


3.03 


" 


2.1 




4.1 


12:06:12.58 


-08 


47:43.12 


2.99 [2.83-3.15] 


2.54 [2.48-2.60] 


2.54 


(2.54) 







4.2 


12:06:11.22 


-08 


47:50.27 


2.54 [2.36-2.71] 


2.35 [2.18-2.64] 


(2.54) 


" 


0.7 




4.3 


12:06:13.15 


-08 


47:59.84 


2.35 [2.16-2.54] 


1.93 [1.91-1.95] 


(2.54) 


" 


0.8 




4.4 


12:06:12.119 


-08 


47:59.51 


-H 


-H 


(2.54) 


" 


1.5 




4.5 


12:06:11.70 


-08 


48:43.04 


3.04 [2.83-3.20] 


2.67 [2.50-3.24] 


(2.54) 


" 


2.3 




5.1 


12:06:13.66 


-08 


48:06.28 


1.87 [1.73-1.99] 


1.76 [1.72-1.89] 


- 


~ 1.85 


0.1 




5.2 


12:06:13.62 


-08 


47:45.85 


1.80 [1.60-1.91] 


1.63 [1.54-1.72] 


- 


" 


1.2 




5.3 


12:06:12.59 


-08 


48:34.06 


1.64 [1.46-1.84] 


1.55 [1.09-1.74] 


- 


" 


0.9 




6.1 


12:06:10.38 


-08 


47:52.10 


2.79 [2.64-2.94] 


2.86 [2.77-2.93] 


- 


~2.8 







6.2 


12:06:10.64 


-08 


48:26.97 


2.65 [2.29-2.79] 


2.53 [2.44-2.71] 


- 


" 


1.5 




6.3 


12:06:11.21 


-08 


47:35.24 


2.77 [2.62-2.92] 


2.53 [2.37-2.64] 


- 


" 


1.1 




7.1 


12:06:15.98 


-08 


48:15.98 


3.55 [3.37-3.91] 


3.47 [3.29-4.04] 


- 


~3.7 





affected by local galaxy weight 


7.2 


12:06:15.95 


-08 


48:17.21 


-H 


-H 


- 


" 





" 


7.3 


12:06:15.95 


-08 


48:18.45 


4.02 [3.70-4.22] 


3.89 [3.34-4.07] 


- 


" 





» 


7.4 


12:06:15.90 


-08 


48:23.01 


-H 


-H 


- 


" 


0.7 


" 


7.5 


12:06:15.75 


-08 


48:27.49 


3.97 [3.78-4.17] 


3.68 [3.43-4.08] 


- 


" 


0.7 


" 


8.1 


12:06:12.33 


-08 


47:28.61 


5.44 [5.19-5.72] 


5.52 [5.15-5.73] 


- 


~5.7 


2.0 




8.2 


12:06:10.58 


-08 


47:49.49 


5.42 [5.17-5.72] 


5.44 [5.10-5.75] 


- 


" 


3.2 




8.3 


12:06:13.25 


-08 


48:03.81 


5.42 [5.17-5.74] 


5.42 [4.91-5.74] 


- 


" 


1.3 




8.4 


12:06:11.36 


-08 


48:44.99 


5.46 [5.21-5.71] 


5.51 [5.22-5.66] 


- 


" 


2.1 


most probable 


9.1 


12:06:11.99 


-08 


47:46.89 


1.74 [1.63-1.85] 


1.74 [1.62-1.78] 


- 


~ 1 


2.7 


candidate system 


9.2 


12:06:11.55 


-08 


47:49.36 


1.76 [1.65-1.87] 


1.75 [1.67-1.79] 


- 


" 


1.7 


" 


9.3 


12:06:12.53 


-08 


48:01.47 


-H 




- 


" 


4.2 


" 


9.4 


12:06:12.31 


-08 


48:01.21 


-H 




- 


" 


1.1 


" 


9.5 


12:06:11.48 


-08 


48:30.61 


1.41 [0.85-1.78] 


1.08 [0.80-1.82] 


- 


" 


5.3 


most probable 


10.1 


12:06:12.13 


-08 


47:44.48 


1.68 [1.48-1.79] 


1.54 [1.50-1.65] 


- 


~ 1.4 


5.7 


candidate system 


10.2 


12:06:11.32 


-08 


47:52.81 


1.57 [1.11-1.72] 


1.60 [1.02-1.76] 


- 


" 


3.5 


" 


10.3 


12:06:12.53 


-08 


48:01.47 


-H 




- 


" 


3.1 


degeneracy with 9.3 


10.4 


12:06:12.31 


-08 


48:01.21 


-H 




- 


" 





degeneracy with 9.4 


10.5 


12:06:11.31 


-08 


48:32.83 


1.38 [0.8-1.86] 


0.94 [0.90-1.07] 


- 


" 


1.6 


most probable 


11.1 


12:06:12.04 


-08 


48:02.25 


-H 




- 


~ 1.3 


6.4 


candidate system 


11.2 


12:06:11.77 


-08 


48:00.49 


-H 




- 


" 


2.9 


" 


11.3 


12:06:13.08 


-08 


48:04.00 


1.42 [1.13-1.78] 


1.72 [1.72-1.72] 


- 


" 


1.9 


" 


11.4 


12:06:12.88 


-08 


47:44.75 


1.12 [0.31-1.48] 


1.10 [0.04-1.15] 


- 


" 





" , bimodal 


11.5 
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Note. — Multipl e-image systems and candidates used and uncovered by our model. Columns are: arc ID; RA and DEC in J2000.0; best 
photo-z using BPZ (Bcnftez 2000 Be nftez et al.||2004| |Coe et al,|2006^ , along with its 95% confidence level; best photo-z using LePhare 
| |Arnouts et al. 1999 llbcrt et al. 2006J, along with its yu% confidence level; spectroscopic redshift, spec-z (images marked in parenthesis 
were not spcctroscopically measured here); z mo deii estimated redshift for the arcs which lack spectroscopy as predicted by the mass model; 
A Position, differenc e in arcsecon ds between the images reproduced by the model and the observed images; comments. System 1 was 
uncovered by |Ebeling et al? ( 2009 ) who measured the redshift of images 1.1 and 1.2 spectroscopically at z = 1.036; o ur spectroscopy yields 
z = 1.033 for this system. Note that we interpret the long arc (images 1.1/1.2) as a double-lensed image, though [Ebeling et al.| | |2009[l 
identify it to consist of several (partially) lensed images. Either interpretation has only a very local (negligible) effect on tne model. All 
other systems listed above are found in our work, where we obtain VLT/VIMOS spectroscopy for systems 1 to 4. Note that unusually 
large errors in the photo-z imply a bimodal distribution, such cases are specified in the comments. We denote in the comments column 
where the most probable image was chosen but other candidates are seen nearby. Systems 9 to 11 are candidate systems, since only some 
of the (< lcr) models can reproduce them, or simply since their photo-z disagrees with the model prediction. Also, systems 9 and 10 show 
similar symmetry to that of system 4, strengthening their identification on one hand, but only 2 radial images are seen. Due to the BCG 
light it is hard to determine unambiguously if, or to which, of these two systems they belong. Note also, that if images 11.1 and 11.2 are 
not multiple images, but, say, a jet coming out of the BCG (see FigurellJ, then images 11.3 to 11.5 may constitute an individual system. 
All other systems we consider as secure identifications, in the context oi the photometric redshifts, internal details, and the reproduction 
by our model. 
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Fig. 2. — VLT/VIMOS spectra obtained for systems 1 to 4, at z sys l = 1.033, z sya 2 y 3 = 3.033, z sys 4 = 2.540. The exposure time for each 
spectrum is shown in parentheses. For the z = 3.03 systems, the pairs are blended into the slits in some cases (2.2+3.2, 2.3+3.3) due to 
less ideal seeing conditions, though in some exposures the seeing was sufficient to resolve each individually. 
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Fig. 3. — Reproduction of systems 1 to 4 by our model, compared 
with the real images. For each system we delens one image to the 
source plane and relens it back to the image-plane to obtain the 
other images of that system. More explicitly, we delens-relens, re- 
spectively, images 1.1, 2.2/3.2, and 4.1, each in the lensing-distance 
ratio expected from its redshift (1.033, 3.03, 2.54, respectively), to 
obtain the results shown here. The upper images in each row are 
the real images while the lower images are those reproduced by our 
model. As can be seen, our model reproduces the images of these 
systems with great accuracy. In these stamp images North is right, 
East is up, and some of the images are slightly zoomed-in for clar- 
ity. Note also, due to their similar position and symmetry, systems 
2 and 3 were considered a single-system in the minimization of the 
mass model. 

distribution and constrain the profile of this cluster, for 
the first time. 

We made use of the location and redshift of all secure 
multiple-images (i.e., excluding candidate systems 9-11) 
to fully constrain the mass model and profile. We ob- 
tained VLT/VIMOS spectra for systems 1 to 4, at red- 
shifts 1.033, 3.03, 3.03, and 2.54, respectively, to help 
pin down the profile with greater accuracy. Details on 
the multiple systems are found in Figure [T] and Table [T 
where in Figure [3] we show examples for reproductions of 
some of the multiple-systems by our model. The % 2 for 
the best model is 25.7. With 35 secure multiply-lensed 
images and 6 free parameters, this yields a reduced x 2 of 
~ 1. Estimating the accuracy of our mass distribution 
quantitatively, for all secure images we obtain an average 
image-plane reproduction uncertainty of 1.3" per image, 
with an image-plane rms of 1.8", typical to paramet- 
ric mass m odels for clusters wit h many multiple images 



(Broadhurst ct al. 2005; Halkola et al. 
|et al.||2007| |Zitrin et al.j|2009b[ ) 
realistic given the ~ 1 noise 



|2006| |Limousin 

Also, this rms value is 

level expected from LSS 



D'Aloisio 



Fig. 4. — Projected total mass profile. The light blue curve shows 
the radial surface mass-density profile in units of the critical surface 
density (re; right-side j/-axis), for a source redshift of z a = 2. The 
dark blue curve shows the overall enclosed mass per radius, M(< 
R) (left-side j/-axis). The widths of the blue curves indicate the 
~ la errors. The thick dash-dotted curve is the (preliminary) 
best-fit model resulting from a 4 X 10 4 step Monte-Carlo Markov 
Chain (MCMC) with Metropolis-Hastings algorithm, which allows 
also the BCG mass to vary. The results are very similar, with only 
some discrepancy at radii below ~ 1". The details of this MCMC 
method will be presented elsewhere. The dash-dotted vertical line 
denotes the mean Einstein r adius (z„ = 2) distribution from 3000 
SDSS clusters at z ~ 0.44 jZitrin et al.||2011b), and the vertical 
dashed line denotes the Einstein radius of MACS1206 (see §4 for 
this comparison). 

of te — 28 ± 3", or ~158 kpc at the redshift of the 
cluster. A projected mass of 1.34 ± 0.15 x 10 14 M Q is 
enclosed by this critical curve for this source redshift 
(see Figure [I]) . For the lower source redshift of system 
1, z s — 1.033, the Einstein radius is ~ 17", enclosing 
a projected mass of 0.8 ± 0.1 x 10 14 M Q . For compari- 
son, our model encloses within 21", a projected m ass of 
M(< 21") ~ l±0.1xl0 14 M Q , while |Ebeling et al.| ( |2009 | 
found similarly, although based on only one system, a 
projected mass of M(< 21") = 1.12 ± 0.05 x 1O 14 M , 
consistent with our result. In addition, as a consistency 
check for the very inner profile around the BCG, we com- 
pared our result to the F160W light. By estimating the 
stellar mass profile from the F160W surface brightness 
photometry, assuming a Kroupa IMF, the stellar mass 
is approximately 10% of the total mass within 30 kpc, 
which is approximately half the BCG effective radius. 
We defer to forthcoming papers the exact assessment of 
the DM mass distribution steepness in the very inner re- 
gion (say, below ~ 5 kpc) after taking into account the 
stellar and gas baryonic contributions. 

The corresponding critical curves for different redshifts 
are plotted on the cluster image in Figure 111 along with 
the multiply-lensed systems. The resulting total mass 
profile is shown in Figure HI for which we measure a slope 
of dlogS/dlogr ~ -0.55 ± 0.1 (in the range [1",53"], 
or 5 < r < 300 kpc; about twice the Einstein radius), 
similar to other usually-relaxed and well - concentrated 



along the line-of-s ight (e.g. Julio et al.| |2010 

& Natarajan 2011 ), the many multiple images used, and 

the small number of free parameters in our modeling. 

For a source at z s = 2.54, the critical curves enclose 
a relatively large area, with an effective Einstein radius MACSi20t> is indeed a relaxed cluster; X-ray and opti 



lensing clust ers (Broadhurst et al. 2005 Zitrin et al.| 



mg 
)b 

' U1 



2009bl 2010). However, it is not perfectly clear whether 
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cal light contours (see e.g., Ebeling et al. 20091 indeed 
show an approximately (circularly) symmetric distribu- 
tion, without a prominent sign of recent merg er. How- 



ever, a high velocity dispersion (1580 km s 1 ; |Gilmour 
et al. 2009), along with the excessive X-ray luminos- 
ity (2.4 x 10 45 erg s " 1 , [0.1-2.4] keV) an d temperature 
(11.6 ± 0.7 keV; see |Ebeling et al.| | 2009b, may imply a 
merg er along the line ot sight (see also Postman et al. 
20111. A full assessment of the degree ot relaxation ot 
this system will be soon enabled, by the dynamical anal- 
ysis from several hundreds member velocities we are cur- 
rently collecting in our spectroscopic program, as well as 
from the combination of other mass diagnostics. 

We note that given its redshift, MACS1206 has a rela- 
tively large Einstein radius. Previous studies have shown 
that other MACS clusters at a redshift of z ~ 0.5 dis- 
tribute around this Einstein radius size [9e — 28"), but 
with a noticeable discrepancy from expect ations of the 
ACDM model and related simu lations (e.g., Zitrin et al. 

even after taking into 



2011a 



Meneghetti et al. 2011) 



account triaxiality-induced lensing bias 

Recently we have applied our lens-modeling technique 
to an unprecedentedly large sample of 10,000 SDSS clus- 
ters, to de duce a representatiy e distribution of Einstein 



radii (see Zitrin et al. 2011b), covering the full clus- 
ter mass-range. Eor the redshift bin corresponding to 
MACS1206 (see Figure 5 therein), the distribution for 



z s = 2 peaks below 10" (with median and mean 9e 
of 4.1"and 5.6", respectively), and rapidly declines to- 
wards larger radii. As may be expected from the dif- 
ferent selection criteria, the Einstein radius of the X-ray 
selected cluster MACS1206 (8 E ~ 26.5" for z s = 2) sits 
at the (far) tail of the distribution: only ~ 1.3% of the 
optically-selected SDSS clusters examined at this redshift 
have similar (or larger) Einstein radii, and as much cor- 
responding enclosed mass, as found in MACS1206. The 
lensing model for MACS1206 presented here, enabled by 
deep, multi-band HST imaging, constitutes an important 
example of the inner mass distributions of systems lying 
at the high-end of the cluster mass function. 
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